Impact of Coulomb interaction and Kondo effect on transport in quantum dots 
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We examine the impact of Coulomb electron-electron interaction on transport in a junction with 
a quantum dot described by Kondo Hamihonian. We analyze the Fermi liquid regime and consider 
the limit of zero temperature. With the help of Keldysh technique we calculate the non-linear cur- 
rent and shot noise as a function of applied voltage. We show that Coulomb interaction markable 
influences the universal effective charge of current-carrying particles e* = 5/3e which can be mea- 
sured in shot-noise experiments. The electron-electron interaction modifies this universal value by 
a factor (e* — 5/3eF) which is less then one and also voltage dependent. 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.63.Kv 



The non-linear conductance and zero-frequency noise 
out of equilibrium (shot noise) can yield information on 
charge fluctuations in physical systems Indeed, in 
many cases the effective charge of carriers may be dif- 
ferent from bare electron charge e. This occurs when 
the fundamental quasiparticle charge does not coincide 
with e or due to interactions. The Kondo effect in the 
quantum dots opens the opportunity to investigate the 
changing of effective charge due to interactions. For this 
purpose it is of a special interest the Fermi liquid regime 
of Kondo model. In this case the backscattering current 
is defined by the scattering term and by the interaction 
quadratic in the spin current. The measurements of the 
backscattering current and shot noise can give the in- 
formation about effective backscattering charge of the 
carriers influenced by interactions. These reasons have 
been put forward in a recent work [l| where the effec- 
tive current-carrying charge of the particles was found to 
acquire the universal value e*/e = 5/3. However, as we 
will show in the present article, the Coulomb interaction 
and corresponding electrical potential fluctuations have a 
noticeable impact on this effective charge. They change 
the value of backscattering charge and cause a voltage 
dependence of e*. To obtain the effective charge with 
or without electrical potential fluctuations we have to 
find both parts (scattering and interaction) of the current 
and shot noise [l|, 0, 0| ■ We concentrate in the following 
on T = limit, where the thermal noise vanishes. We 
consider the regime around the strongly interacting fix 
point (Fermi liquid regime) of Kondo Hamiltonian, that 
is, eV << Tk {Tk is the Kondo temperature) counting 
for potential fluctuation at the one loop approximation. 
Here we notice that the role of Coulomb interaction in 
the case of a large dots was analyzed before in a number 
of articles 0, Q ■ 

Our principal result can be written compactly as 
e*/e = |F, where factor F depends on the bias voltage, 
charging energy Ec = e^/C (here C is the total capaci- 
tance), and the effective number of conducting channels 
N . In Fig.l we plot this factor as a function of voltage 
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FIG. 1: Voltage dependence (in units of Coulomb energy Ec) 
of factor F. The curves are related to different number of effec- 
tive conducting channels N starting from N=2 up to N=1Q 



for different N . For validity of perturbation theory that 
we use, the absolute value of the voltage must be not too 
small. Also for dot systems Eq » Tk- In this limit we 
have 
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In the region of small voltages factor F (Figl.) rather 
strongly ( 30%) reduces the effective charge. Now we 
present the basic steps which lead to the derivation of 
equation ([T]). 

Effective action: Coulomb interaction in a quantum 
dot is responsible for formation of Kondo effect and for 
electric potential fluctuations in the whole junction. The 
Kondo effect depends also on the tunnelling rates be- 
tween the dot and the leads. The later strongly influence 
the Kondo temperature Tk- The Coulomb interaction is 
characterized by the charging energy Ec- The parameter 
that controls the impact of electron electron interactions 
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on the tunnelling is the total conductance g or the num- 
ber of open channels N . More precisely, in the case of 
large N we can find the interaction corrections by per- 
forming 1/iV expansion. Here we address the unitary 
limit of Kondo quantum dot with odd numbers of elec- 
trons in it, that is, we deal with ordinary Kondo effect 
in a quantum dot. Though, in this case N = 2 (effec- 
tive N may be increased by a shunting resistance) we 
believe that even at 1/A^ level the qualitative picture is 
correct. The problem is also interesting because of its 
relation to a rare case of a non-equilibrium system with 
interaction where the electrical potential fluctuations are 
included. As an advantage of the unitary limit is that 
the effective transmissions are close to the unity. Indeed, 
conductance in additions to (Jq = 2e^//i has a small non- 
linear contribution ay ~ —V'^/T'^, where the bias volt- 
age V « Tk- As we will show the Coulomb interaction 
modifies only the nonlinear part of the transport current 
and shot noise. Therefore, to the lowest order in V/Tk 
one can calculate the effective action for electric poten- 
tials taking in considerations only the unitary limit. In 
this limit we can use two approaches which give the same 
result for spin 1/2 Kondo model: one is based on scat- 
tering representation of the Kondo Hamiltonian Q, the 
other is the mean field slave boson approximation (SB) 
0. The SB is more transparent and it will be used to 
obtain the effective action. We start with the Anderson 
Hamiltonian 

H = ffi + fffi + ^(i;„cL,A + H.c.) 

+ Y,edld, + Ud\d^d\d^ (4) 

a 

The first two terms correspond to non-interacting elec- 
trons in the two leads 
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where Caa;k, ^ak are the electron field operator and the 
electron energy of a lead. Index a = L,R indicates left 
(right) lead. We assume that the leads are biased by 
the fluctuating electric field potential drop v{t) = . 
There are also electron electron interactions in the leads 
which are not included in the Hubbard coupling U. At 
first we separate the Hubbard interaction that leads to 
Kondo effect. For this we apply to Hamiltonian ^ the 
SB approximation at f7 — ^ oo. In the slave-boson ap- 
proach, the localized electron operator dj. is represented 
by fib with b and /J being, respectively, the standard bo- 
son and fermion operators. The action on Keldysh con- 
tour then will incorporate the Anderson Hamiltonian (|4]) 
with the new boson and fermion operators. The quadric 
term is replaced by the two new terms: one represents 
the renormalized level position in the dot e, and the sec- 
ond stands for requirement of a single occupancy which 



demands a Lagrange multiplier A 

b^b + J2PaU = '^ 



(6) 



In the SB approximation the Bose operators 6^ , b are re- 
placed by their expectation value b. Having the Anderson 
part of the action arranged we can get the general form 
of the total effective action by including the Coulomb 
interaction term. Performing the standard Hubbard- 
Stratonovich decoupling of this term in the action, one 
introduces on the Keldysh contour fluctuating electrical 
potentials Vi{t) = ^^,(i = 1,2) Next we in- 

tegrate out the fermion operators in the leads, and also 
add the action of external circuit expressed in terms of 
shunt resistance Rg- Thus the effective action acquires 
a form Sef / — Sb + S + Ss, where nonoperator bosonic 
part of the action is 



Kib^~l)a: 
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where CTz is the Pauli matrix. The action S consists of two 
parts: the one describing the scattering of electron in the 
leads, and the other {Shs) responsible for the fluctuating 
Habbard-Stratonovich fields 



5* = -iTrLnG- 



Shs 



(8) 



G-\t,t') - [^--i-—a,{h + ~9R)]<y. (9) 



~9LMt,t') = e±'^«5i,fl(i-<')e^^^^''\ 

C f ■ ■ 

Shs = TT^ / dt(tripazip) 
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dtdt'xithit - t')x{t') 



+ J dtx{t){eV~{^i+^2)/2)] (12) 

here F is a tunnelling width, x = V^i — ^2, the Fourier 
transform of ^{t) at finite temperature T is 7(w) = 
Locoth{uo /2T). Here and below the values with the hat 
w, mean the diagonal 2x2 matrices with entries wi, V2 
and (^1, </32, respectively. gL,R represent Keldysh matrices 
of the electron propagators in the leads at the position of 
the dot. Their Fourier transforms are known: = —in, 
(e) = 2i7r/(e) with /(e) as Fermi distribution func- 
tion. All dependencies on constant bias and fiuctuating 
potentials are collected in the exponents of Eq. pU]) . We 
also define Tk — Tb^ as an effective Kondo temperature. 
The G{u}) includes the Lagrange multiplier which shifts 
the localized level position e to e = e-|-A. For our purpose 
it is enough to find both free parameters, and renor- 
malized level e, by solving two self-consistent equations 
in equilibrium (neglecting the fluctuations of electrical 
potential) [9]. These equations follow from extremum 
conditions of effective action Sef / relative to b and e 

dSeff _ dSeff 



db 



de 







(13) 



3 



Unlike the multichannel Kondo model, the solution 
for spin 1/2 Kondo quantum dot satisfies @ e ~ 
Tk{Tk/^) « Tk which permits us to use e/Tx as a 
small parameter. 

To obtain the effective action at one loop level we have 
to expand the logarithm in equation ([5]) to the second 
power in quantum field x- The calculations are straight- 
forward, though, rather long. The results for the unitary 
limit are in accordance with those obtained by bosoniza- 
tion technique 0, |Tl|, [13, [3] and in reference [3| • The 
part of effective action (~ Ln) ([5]) which describes the 
electron scattering thus acquires a form 



dtdt'x{t)^{t-t')x{t') 



dtx{t){^i+^2)/2] 



(14) 



Here the number of the channels N = 2. Collecting to- 
gether all terms we arrive to a compact form of effective 
action which is the principal result of this section: 



J J dtdt'ip{t)B-\t-t')ip{t') 
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The Fourier image of Coulomb Green function B(t) we 
find B(u;) = B^(w)/a;^ where 



(16) 



here Q = EcN/n, N = N{1 + h/{Ae^Rs)) is the effective 
number of conducting channels, and 7 = 7(w). Bt,(a;) 
is matrix Green function (GF) of electric potentials Vi. 
Having the effective action we can calculate correlation 
functions of fluctuating potential flelds and, thus, find 
the conductance and shot noise corrections in the Fermi 
liquid regime of the Kondo model. 

Fermi Liquid Regime: Near the unitary limit it is con- 
venient to follow the scattering approach. For this pur- 
pose we use the basis of s and p scattering states rather 
than those of the left-lead and right-lead states. At 
V = Q the p states are decoupled from the dot. For non- 
zero fluctuating potentials the p-states are not decoupled, 
however, the s, p basis is important. In real junctions 
with strong Kondo effect, nearly T = the dot is de- 
scribed by strongly interacting flx-point with many body 
state (5). Before scattering the phases of both states, h 
and s -state, coincide, while passing the scattering region 
h state acquires only an extra phase tt compared to that 
of s-state. Then the fix-point Hamiltonian can be written 
[1, H, in the new basis {ba,k,a,a.k) (here a^.fc stands 
for p-states) as H = Hq + Hg + Hint 

+ l^iblkG^^k + alj^b^^k) 

k.cr 



H, 



{a + ^k')bl,b,^k' (17) 
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where b = 2a and a — l/(27r) [14] and — J2k ^T-fe- The 
Kondo temperature Tk is the only energy scale of the 
fixed-point Hamiltonian. We can integrate out the Oo-./t- 
states and obtain the effective action on Keldysh contour 
only for interacting states b 



S = 



dtdt'J2bl{t)G;,'it 

ka 



t')bka{t') 



dt{Hl + Hl^^a: (19) 

Glkii^^') = [[i^^-ik)<y.-\v<y.Gak<yzv] (20) 

Here Gak is matrix Green function for decoupled nonin- 
teracting a- states (for example, G^^,(e) = l/(e — ^fc-f i(5)). 
We separate potentials in (^0)1 into fiuctuating part and 
constant applied voltage F , that is, write v ^ V + v. To 
obtain the Coulomb interaction corrections to conduc- 
tance and shot noise we have to express the current in 
terms of Green functions (the inverse of ([20|l and then, 
with help of effective action , to average the product 
of fluctuating potentials. This goal can be achieved by 
applying the perturbation theory (it corresponds to 1 /N 
expansion) . 

Current: To begin with, we consider, at first, the aver- 
age current. The backscattering current operator is given 
as daii: 
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1^2vTk 
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and ricr = 6^60-, ct = — cr. The sum over a stands for spin 
summation, which in our case is trivial. The averaged 
backscattering current consists of two parts h = Is+hnt, 
where Is and lint originate, respectively, from the scat- 
tering Hamiltonian Hs and the interaction Hint- Is and 
lint are generally expressed in terms of Green's functions 
for a and b states. 

= -|(^)' E(^fc+^'^')'(^fe,fe'-^i2^) 
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= < [GakVCTzGbk^y zGbkiYtt > 
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> 



< [Gbki<yzGbkVcrzGak 

< [GakV<JzGbk<JzGb]u[GbGb]u > 

< [Gb<y zGbkV(TzGak]u[GbGb]\l > 



(23) 

(24) 
(25) 
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FIG. 2: Diagrams for processes which contribute to Is- The 
solid lines with momentum k denotes 6-electron Green func- 
tions gbk{e). Dashed lines represent GF of noninteracting 
states a. The dot line stands for Coulomb GF iBv(w) (|16|) . 
All Green functions are 2x2 Keldysh GF. Bias voltage V 
denotes the vertex that connects dash and solid line. Each 
vertex (but most left one) includes Uz, and matrix product is 
considered. 



where brackets < ... > denote the average over fluctu- 
ating potentials with the effective action The ex- 
pression in square brackets is matrix product in Keldysh 
space and convolution in time. One has to take the 11 
component of this matrix product. The last expression 
(|23)) . unhke the scattering part of the current ((22)) . in- 
cludes, in addition, the Green functions at the dot po- 
sition Gb = '^j^Gbk- In the zero order approximation 
(no potential fluctuations) for the Fermi liquid regime 
of Kondo effect the nonlinear current was found before 
(see also recent works [H, 01 j where the ex- 
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phcit separation between scattering and interacting con- 
tributions to conductance was done). In this case the 
backscattering current reads as. 



3h 



(26) 



Let us now include the fluctuations of potential flelds. At 
first we do the calculations for Is. For this purpose we 
expand Gb to the second order in fluctuating electrical 
potentials. The GFs for 6-states to zero approximation 
in fluctuating flelds acquire simple form: 



7^^ 
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V/2 + iS Lu ~^k + V/2 + i6' 
(1-/(&))(5''(M-5^(M) 



9bk{^) 
9ll{^) 
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The diagrams are a convenient way to represent different 
contributions to the transport current. For Ig two types 
of diagram are shown in Fig2. The principal contribu- 
tions to the scattering part of the current at 1 /N order 
are given by a number of processes. These processes are 
presented by: (i) two diagrams in Fig2, (ii) two diagrams 
similar to those in Fig2, however, with the Coulomb (dot) 
line connected with other vertex (V) (these diagrams are 
obtained by an interchange between the V vertex and the 




FIG. 3: The processes that contribute to the interacting 
part of transport current Unt- The notations for momentum 
dependent GF are the same as in Figl. The lines without 
momentum denote integrated over momentum GFs gb- 



end of the dot fine ), and, (ui) four "conjugate" diagrams 
for Ak^k' (like the cases (i) and (ii)). These "conjugate" 
diagrams together with diagrams in Fig2 simply double 
the real parts of those in Fig2. We notice that there are 
yet processes which can be find by self-energy insertion in 
the one momentum line. However, their impact on cur- 
rent is small. Thus summing all relevant contributions 
and including the part without fluctuations one finds the 
scattering part of the total current 



9 2Q 
^ In 

16iV^ eV^ 



+ 1)] (27) 



We need more diagrams (see Fig3) to calculate lint- Di- 
agrams of types (a)and (c) are natural generalization of 
those in Fig2. Diagrams of type (b) and (d) are of "scat- 
tering" type and appear only due to the interactions. 
The accounting for relevant processes goes as follow: at 
first the diagrams for fl are "conjugate to those of II and 
after summation, as in the case for 1^, they double the 
real parts of diagrams in Fig3. Next, we have :(i) four 
diagram of type (a) (two of them appear due to inter- 
change with vertex V) , (ii) 2 4 = 8 diagram similar 
to (c) (factor four appears due to four ways to connect 
V vertices), (iii) two diagrams of type (b), and (iv) four 
diagrams of type (d) . The summation ah of this contri- 
butions results in the Coulomb interaction correction to 
the lint 

/». = ^^4(^f[l + ^(ln[^]+0.11)] (28) 

The corresponding corrections to differential conduc- 
tance due to electron-electron interaction are obtained 
by taking the derivative of (I27p . ([28|) with respect to V. 

Shot Noise: The effective charge of the current- 
carrying particles is deflned as e* = Sn/hi where Sfi is 
the zero frequency Fourier transform of current-current 
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which includes potential fluctuations, we obtain 




FIG. 4: Diagrams which contribute into the scattering part 
of the noise. Notations of the GFs and momentums are the 
same as in Fig2, however, the time contour indices of the edge 
vertices are different: here index 2 is in the left and index 1 
in the right vertex. Diagrams b) enter with sign minus. 




FIG. 5: The processes that contribute to the interacting part 
of shot noise The notations for GF correspond to those 

in Fig2 and Fig3. The edge Keldysh vertex indices are the 
same as in Fig3. Diagrams b) enter witli minus sign. 



correlation function, 



S{t) =< I{t)I{0) > - < / >2 



(29) 



The same as the current, noise power can be separated 
into the two terms S = + 5*™* related, respectively, to 
the scattering part of the Hamiltonian Hg and to the in- 
teraction part Hint ■ In the Fermi liquid regime (without 
fluctuating potentials) at temperature T = the shot 
noise was obtained before [J, Q 



h ^2 



3 > Tl 



(30) 



The influence of Coulomb interaction changes the scat- 
tering and interacting parts of Sn- Explicitly, we can 
find these changes (like in the case of the current) at 
order of perturbation theory. In the formula for 
noise ([50)) we assign Keldysh indices to the current oper- 
ators pT|) . and then, using Keldysh technique, sum the 
all relevant diagrams. For scattering part these dia- 
grams are shown in Fig. 4. There are also equivalent (not 
represented by Fig. 4) diagrams. They also contribute to 
the noise, and their number, like for current, can be eas- 
ily counted. Thus, for the total scattering noise power 



i(ln[^] + 1.18)] (31) 



For interacting part of the shot noise 5'™* the typi- 
cal diagrams are represented in Fig. 5. In addition to 
them there are also "scattering" diagrams (GFs of inside 
bobble point in one direction). Each such diagram cor- 
responds (with the same contribution) to one in Fig. 5. 
After summation of all diagrams we obtain 



Qtnt 



2h 



(^)2[H--^(ln[^] 



-0.31)] (32) 



By using Eqs. pil32l27l28p and definition of the e* we 
can immediately obtain factor F (see EqllJ. 

Conclusions: We have presented a theory to evaluate 
the effects of electron-electron interactions on quantum 
transport through the Kondo quantum dot. We address 
the important problem of effective backscattering charge 
for current-carrying particles in the interacting systems. 
This charge can be obtained from the measurements of 
shot noise and backscattering current in the Fermi liquid 
regime of a Kondo quantum dot. We show that in this 
regime the potential fluctuations can markable reduce 
the recently calculated universal value 5/3e 0], though, 
the effective charge remains bigger then e. To get our 
principal result we have calculated the transport current 
and obtained the zero frequency shot noise power. Of 
cause, the potential fluctuations are not the only cause 
to change the effective charge. The pure unitary limit 
in real system may by broken by the symmetry break- 
ing interactions like potential scattering or magnetic field 
which can also influence the effective charge. However, 
in the case these effects are small, they can be properly 
estimated P, i Q [H and render a small correction to 
e*. 
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